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Abstract 



We consider the possibility of observing CP violation in quasi-inclusive decays of 
the type B~ K~X, B~ -> K*~X, B° K~X and B° -> K*~X, where X 
does not contain strange quarks. We present estimates of rates and asymmetries for 
these decays in the Standard Model and comment on the experimental feasibility 
of observing CP violation in these decays at future B factories. We find the rate 
asymmetries can be quite sizeable. Observation of such asymmetries could be used 
to rule out the superweak model of CP violation. 
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I. INTRODUCTION 



The possibility of observing large CP violating asymmetries in the decay of B mesons motivates 
the construction of high luminosity B factories at several of the world's high energy physics labora- 
tories. The theoretical and the experimental signatures of these asymmetries have been extensively 
discussed elsewhere 0, @, @, 0, ||. At asymmetric B factories, it is possible to measure the 
time dependence of B decays and therefore time dependent rate asymmetries of neutral B decays 
due to B — B mixing. The measurement of time dependent asymmetries in the exclusive modes 
B° — > ipK s and B° — > 7r + 7r~ will allow the determination of the angles in the Cabbibo-Kobayashi- 
Maskawa (CKM) unitarity triangle. This type of CP violation has been studied extensively in the 
literature. 

Another type of CP violation also exists in B decays, direct CP violation in the B decay 
amplitudes. This type of CP violation in B decays has also been discussed by several authors 
although not as extensively. For charged B decays calculation of the magnitudes of the effects 
for some exclusive modes and inclusive modes have been carried out ||, [^], ||, |§, HlO] , |T 



12| . In contrast with asymmetries induced by B — B mixing, the magnitudes have large hadronic 
uncertainties, especially for the exclusive modes. Observation of these asymmetries can be used to 
rule out the superweak class of models [[OJ . 

In this paper we describe several quasi-inclusive experimental signatures which could provide 
useful information on direct CP violation at the high luminosity facilities of the future. One of the 
goals is to increase the number of events available at experiments for observing a CP asymmetry. In 
particular we examine the inclusive decay of the neutral and the charged B to either a charged K or 
a charged K* meson. By applying the appropriate cut on the kaon (or K*) energy one can isolate a 
signal with little background from b — > c transitions. Furthermore, these quasi-inclusive modes are 
expected to have less hadronic uncertainty than the exclusive modes, would have larger branching 
ratios and, compared to the purely inclusive modes they may have larger CP asymmetries. In this 
paper we will consider modes of the type B — > K(K*)X that have the strange quark only in the 
K(K*)-meson. These processes include contributions from the one loop process b — > sg* — > sqq as 
well as the tree level process b — > uus. The interference between these two processes is responsible 
for the direct CP violation. 

In the next section, we describe the experimental signature and method. We then calculate 
the rates and asymmetries for inclusive B~ — > K~(K*~) and B° — > K~(K*~) decays. In the last 
section, the theoretical uncertainties in the calculation are discussed. 



II. EXPERIMENTAL SIGNATURES FOR QUASI-INCLUSIVE B -> SG* 

In the T(4S') center of mass frame, the momentum of the from quasi-two body B decays 

such as B -> may have momenta above the kinematic limit for mesons from b — ► c 

transitions. This provides an experimental signature for b — ► sg*, g* — > uu or g* — > dd decays where 
g* denotes a gluon. This kinematic separation between b — > c and b — > sg* transitions is illustrated 
by a generator level Monte Carlo simulation in Figure 1 for the case of B — > K*~ . (The B — > K~ 
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spectrum will be similiar). This experimental signature can be applied to the asymmetric energy 
B factories if one boosts backwards along the z axis into the T(4S) center of mass frame. 

Since there is a large background ("continuum") from the non-resonant processes e + e _ — » qq 
where q = u,d, s, c, experimental cuts on the event shape are also imposed. To provide additional 
continuum suppression, the ll B reconstruction" technique has been employed. The requirement that 
the kaon and n other pions form a system consistent in beam constrained mass and energy with a 
B meson dramatically reduces the background. After these requirements are imposed, one searches 
for an excess in the kaon momentum spectrum above the h — > c region. Only one combination per 
event is chosen. No effort is made to unfold the feed-across between submodes with different values 
of n. 
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FIG. 1. Generated inclusive B — * K*~ momentum spectrum. The component below 2.0 GeV/c is due 
to b — > c decays while the component above 2.0 GeV/c arises from quasi-two body b — > sg* decay. The 
normalization of the b —* c component is reduced by a factor of approximately 100 so that both components 
are visible. 



Methods similar to these have been successfully used by the CLEO II experiment to isolate 
a signal in the inclusive single photon energy spectrum and measure the branching fraction for 
inclusive b — > 37 transitions and to set upper limits on b — > scf) transitions |]14 |, | 15| . It is clear 
from these studies that the B reconstruction method provides adequate continuum background 
suppression. 

The decay modes that will be used here are listed below: 

1. B~ -»■ K^-tt° 

2. B° -> K^-n+ 
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3. 


B~ - 


-»■ K(*> 


7T 


7T~ 


h 


4. 


B°- 


-> 


7T 4 


-7T° 




5. 


B° - 


■* tfW" 


7T 4 


~7T~ 


7T+ 


6. 


B~ - 


-> tfW- 


~7T" 


f 7T- 


-7T° 


7. 


B~ - 


-»■ K { *> 


~7T" 


f 7T- 


~7T + 7T~ 


8. 


B°- 


-> #w- 


7T 4 


~7T~ 


7T + 7T 



In case of multiple entries for a decay mode, we choose the best entry on the basis of a x 2 formed 
from the beam constrained mass and energy difference [i.e. \ 2 = (Mb/SMb) 2 + (AE /5AE) 2 ]. In 
case of multiple decay modes per event, the best decay mode candidate is picked on the basis of 
the same \ 2 - 

Cross-feed between different h — > sg decay modes (i.e. the misclassification of decay modes) 
provided the is correctly identified, is not a concern as the goal is to extract an inclusive 

signal (an example of the signal is shown in Fig. |2|). The purpose of the B reconstruction method 
is to reduce continuum background. As the multiplicity of the decay mode increases, however, the 
probability of misrecontruction will increase. 

The signal is isolated as excess K^*'" production in the high momentum signal region (2.0 < 
Pk(*) < 2.7 GeV) above the continuum background. To reduce contamination from high momentum 
B — > 7r~(p~) production and residual b — > c background, we assume the presence of a high mo- 
mentum particle identification system as will be employed in the BABAR, BELLE, and CLEO III 
experiments. 

We propose to measure the asymmetry iV(if ( * )+ - K ( -*^~)/N(K^ + + K^~) where if(*) ± orig- 
inates from a partially reconstructed B decay such as B — » (mr)° with 2.7 > p(.?r < -*- ) ~) > 2.0 
GeV. The additional pions have net charge 0, n < 4 and at most one neutral pion is allowed. We 
assume that the contribution from B K~rj X decays has been removed by cutting on the rj 



region in X mass. It is possible that the anomalously large rate from this source |16j could dilute 
the asymmetry. 



III. EFFECTIVE HAMILTONIAN 

In the Standard Model (SM) the amplitudes for hadronic B decays of the type b — > qff are 
generated by the following effective Hamiltonian ||T7|: 



Hi ff = ^[v fb v f * q ( Cl oi f + c 2 oi f ) - + v cb v; qC c + V tb V*cl)Ol] + H.C. , (1) 

y 1 i=3 

where the superscript u, c, t indicates the internal quark, / can be u or c quark, q can be either a 
d or a s quark depending on whether the decay is a AS = or AS" = — 1 process. The operators 
Of are defined as 
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FIG. 2. Monte Carlo simulation of the inclusive B — * K*~X signal with the B reconstruction method: 
(a) the beam constrained mass distribution (b) the distribution of energy difference (c) the K~7r° invariant 
mass after selecting on energy difference and beam constrained mass 



O q fl = q a l»Lfpf^Lb a , 0\j = ^Lff-fLb , 

3 % = qi,Lbq'rL(R)q' , Of, 6 = q al ,Lb^L(R)q' a , (2) 
01, = ^Lbe q ,qYR(L)q' , Of, 10 = ^-q al ,Lb p e ql q'^R{L)q' a , 

where R(L) = 1 ± 75, and q' is summed over u, d, s, c and b. O2 and 0\ are the tree level and 
QCD corrected operators, respectively. O3-Q are the strong gluon induced penguin operators, and 
operators O 7 _i are due to 7 and Z exchange (electroweak penguins), and "box" diagrams at loop 
level. The Wilson coefficients c{ are defined at the scale /i ~ m;, and have been evaluated to next- 
to-leading order in QCD. The c\ are the regularization scheme independent values obtained in Ref. 
||. We give the non-zero c{ below for m t = 176 GeV, a s (mz) = 0.117, and /1 = = 5 GeV, 

Cl = -0.307 , c 2 = 1.147 , 4 = 0.017 , c* = -0.037 , c* = 0.010 , c* = -0.045 , 
4 = -1.24 x 10" 5 , 4 = 3 - 77 x 10_4 > 4 = -0.010 , 4o = 2- 06 x 10 ~ 3 , 
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cZ = -cZ/N c = Pr/N Cj c u 7 % = Pr, cZo = (3) 

where N c is the number of color. The leading contributions to P l se are given by: P l s = (f 2 -)ca(y + 
G(m h n,q 2 )) and P\ = (^ L )(N c c 1 + c 2 )(^ + G{m h ^ ) q 2 )). The function G(m,[i,q 2 ) is given by 

G(m,/i,g 2 ) = 4 / x(l - x)ln ^ — dx . (4) 

All the above coefficients are obtained up to one loop order in electroweak interactions. The 
momentum q is the momentum carried by the virtual gluon in the penguin diagram. When q 2 )4m 2 , 
G(m, n, q 2 ) becomes imaginary. In our calculation, we use m u = 5 MeV, = 7 MeV, m s = 200 



MeV, m r = 1.35 GeV 18.19 



We assume that the final state phases calculated at the quark level will be a good approximation 
to the sizes and the signs of the final state interaction (FSI) phases at the hadronic level for quasi- 
inclusive decays when the final state particles are quite energetic as is the case for the B decays in 
the kinematic range of experimental interest || . As pointed out by Gerard and Hou |7[] and clarified 
by Wolfenstein ||. when calculating rate asymmetries using the absorptive amplitude given above, 
one must be careful to be consistent with the requirements of the CPT theorem. Gerard and 
Hou |7j noted that CPT is violated if one does not include all diagrams of the same order. The 
interference term responsible for the rate asymmetry due to Cj w contains two contributions: an 
interference between penguin amplitudes of order a 2 and a contribution from the interference of 
the tree amplitude with a higher order penguin diagram that contains an absorptive part from a 
vacuum polarization bubble in the gluon propagator. These two contributions cancel each other. 
Therefore, in practical calculations Cj" must be treated as real. The general rule is that the phase 
of the penguin Wilson coefficient must be dropped if there is a tree amplitude with the same CKM 
factor and the final states for the tree and penguin amplitudes are the same. In a more general 
analysis of this problem from CPT and unitarity consideration ||, Wolfenstein showed that diagonal 
strong phases (the phases due to the rescattering of the state which is the same as the final state 
e.g. uu — > uu) do not contribute to partial rate asymmetries. The phase in Cj n is a diagonal phase 
in this sense. We will follow this prescription of Ref. M to remove the redundant strong phases. 



IV. MATRIX ELEMENTS FOR B~ -> K~X AND B° -» K~X 

We proceed to calculate the matrix elements of the form (KX\H e jf\B) which represents the 
process B — > KX and where H e ff has been described above. The effective Hamiltonian consists of 
operators with a current x current structure. Pairs of such operators can be expressed in terms of 
color singlet and color octet structures which lead to color singlet and color octet matrix elements. 
In the factorization approximation, one separates out the currents in the operators by inserting the 
vacuum state and neglecting any QCD interactions between the two currents. The basis for this 
approximation is that, if the quark pair created by one of the currents carries large energy then it 
will not have significant QCD interactions. In this approximation the color octet matrix element 
does not contribute because it cannot be expressed in a factorizable color singlet form. The color 
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octet operators could contribute if, for instance, the quark pair emits or absorbs a gluon [ 2IJ . It has 



been shown |2TJ that in the leading order, where the energy of the light quark pair E ~ with 
mb — > oo, the octet matrix element vanishes when the final state is two 0~ mesons. In our case, since 
the energy of the quark pairs that either creates the K or the X state is rather large, factorization is 
likely to be a good first approximation. To accommodate some deviation from this approximation 
we treat N c , the number of colors that enter in the calculation of the matrix elements, as a free 
parameter. In our calculation we will see how our results vary with different choices of N c . The 
value of N c ~ 2 is suggested by experimental data on low multiplicity hadronic B decays ||. 

In the factorization approximation the matrix element of B~ — > K~X decay can be expressed 

as, 

M = M 1 + M 2 (5) 



where 



where 



Ml = ^(K-\sY(l-l 5 )b\B) £ (X|g 7 M{^(l-7 5 )+^(l + 7 5 )}?|0) 



M 2 = ^(FL U {X sy (1 - 7 5 ) b B)(K~\ sy»(l - 7 5 ) u | 0) 
v2 

+FB„(A' | u (1 - 7 5 ) b B)(K-\s(l + 7 5 ) u | 0)) (6) 



L « = v "( Cl + w) +A > + W A < + A ' + W Am 

L ' =A " + w c Ai -\i A>+ k Ai0 ) 

L, = A 3 + j-A,- 1 -(A, + l-A w 

R « = A S + W A <> + A T + -k A * 

R ' =A " + w c At -\i A7+ k A> 
R ° =A ' + k Ai -\( A7+ k A> 

FR„=-2 (^As + Ae + ^Aj + Az) 17) 



We have defined 



q=u,c,t 
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with 



v q = v q * s v qb 



Using the definition 

(K-\sY(l- 1 5 )u\0) = if KP » K 
where fx is the kaon decay constant, and the free quark equation of motion one has 

(K-\s(l + 1 5 )u\0) = - i - fKm2K 



m u + m s 

Using these two results we can simplify M 2 and write it in the form 

M 2 =if K [a(X\u{l + ^)b\B-)+ (3(X\u(l--f 5 )b\B-)] 

with 

a = m b FL u 

FR u m\ 

13 = —m u FL u ■ 

m s + m u 

To calculate Mi we express 

(K-\ 7sY(i -i 5 )b\B) = u( PB + PK y + /_(p B - PK y 

where /+, /_ are Lorentz invariant form factors which are functions of (ps — Pk) 2 ■ 
For the decay B° — > K~X decay, M\ = and only M 2 contributes. 

V. MATRIX ELEMENTS FOR B~ -► K*~X AND B° -► K*~ X 

For B~ — > K*~X decay, we also write the matrix element as 

M M • M> 

where Mi has a similar structure as in Eq. (6). M 2 has the form 

M 2 = FL u (X\uY(l - 7 5 )&|£>™*-0*.eJ 

where 

(K x *\s 1 "(l- 1 5 )u\0) = m K .g K .e ll x * 



with gx*, the decay constant and e A1 A * being the polarization vector of the vector meson. 

For the B° — > K*~X decay only M 2 contributes. To calculate Mi, following the notation of 
Bauer, Stech and Wirbel (BSW) p2j we write 

(K*~\ s-y»(l - 7 5 ) 6 | B-> = he^e^pl* (18) 

+i{6 2 £ M * - & 3 (pb + • k + & 4^* • kk„} (19) 



with 



= -^*L (20) 

6 2 = (m B + m K ,)A l {k 2 ) (21) 

b 3 = (22) 

h = 2m K ,[A (e) - A 3 (fc 2 )] 



where k = ps — Pk*- In our calculation we will use the form factors of Ref p3 for the primary 



result. To check the dependence of the results on form factors we will also use the modified BSW 



model |24] which has a dipole behavior for the form factors. The form factors in Ref [p3| are first 
evaluated at k 2 = and then extrapolated to a finite k 2 assuming a monopole behavior for all the 
form factors. 

We considered the possible contribution from annihilation graphs to both decay rates and asym- 
metries. In agreement with previous estimates, the annihilation contribution to rates is found to 
be small [|25[| . The contribution to CP asymmetries is potentially interesting since the dependence 
on CKM parameters is quite different. In this case if we limit ourselves to the processes b — > suu, 
b — > sdd that have only a strange quark in the meson then the contribution to the asymmetry 

from the annihilation term turns out to be too small to be of interest. 



VI. DECAY DISTRIBUTION AND CP ASYMMETRIES 

In this section we describe the formalism to calculate the decay distribution, asymmetries and 
the decay rates. 

The general form of the matrix element is 

M = Mi + M 2 (24) 

and so 

|M| 2 = (Mxl 2 + |M 2 | 2 + M\M 2 + M X M\ (25) 
Now I Mi 1 2 has the structure 
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where 



and 



M 1 \ 2 = H /X1/ W>" (26) 



H, v = (K-, (K*~)\ J, \B-)(B~\ Jj \K~, (K*~)) (27) 



= J2(^)^(Pb-Pk-Px)(0\J,\X}(X\JJ\0) (28) 



with 



J, = E ?7m{A,(1 - 7 5 ) + ^(1 + 7 5 )}g (29) 

u,d,s 

In the parton model approximation we can interpret the above process as the decay 

B(p B ) - K(p K ) + q( Pl ) + q(p 2 ) (30) 



with p x = Pi +P2- 
We can also express 



= 2ImtJ d^xe-^iO \T[J^x) jt(0)]| 0) (31) 



with q = Pb ~ Pk — Px- The parton model approximation is the leading term in the expansion for 
the T product in the above equation and so this form for is useful to calculate higher order 
corrections to the parton model approximation. 

The decay distribution is given by 



dT 1 1 



Mtfdmiz (32) 



dE K (2tt) 3 16m| 
where m\ 2 = (p 1 + p K ) 2 and \Mi\ 2 has the structure 

|M X | 2 = 24 E {[(l^l 2 + \R q ?)K» + (\L q ? ~ \R q \ 2 )Al]H^ - 2m 2 q g, v H^Re{L q R* q )} (33) 



with 



A\ v = Pl^v + P\vP2p - g^vPl ■ P2 

A% = te^p 2a p^ (34) 

(35) 



(See the Appendix for the full form of |Mi| 2 .) 
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For decays involving K , |M 2 | 2 has the form 

|M 2 | 2 = | (X I u {a(l + 7 5 ) + P(l -l 5 )}b\B) \* (2tt) 4 S 4 (p B - p K - p x ) (36) 

x 

= ^(2tt) 4 5 4 ( Pb - p K - p x ) (B | J | X) (X | Jt | B) 



with 



J f =u{a(l + 7 5 )+/3(l- 7 5 )}6 



In the parton model approximation we replace 

d 3 p 



EWhE/ j^^H Pu ,s))(u( Pu ,s)\ 



(37) 



(38) 



(39) 



where \u(p u , s)) is a free quark state with momentum p u and spin s. As in the previous case it is 
also possible to express \M 2 \ 2 as 



|M 2 | 2 = 2Im{B\i J d 4 xe ipKX T[J{x)j\0)]\B) 



(40) 



where the parton model approximation is again the leading term in the expansion of the T product 
in the above equation and can be interpreted as the two body process b — > Ku. 

In the parton model approximation we can write for K~ decay, 



|M 2 | 2 = 4/ 2 (|a| 



l )Pb ■ Pu + 2Re(a(3*)m b m u 



and for K* decay 



I Mo 



j I 7-1 |2 2 2 



Pb'Pu + 



2pb ■ Pk*Pk* ■ Pu 



(41) 



(42) 



For the interference terms, we have, for the K decay 

M X M\ = -if K (K-\sY(l - l b )b\B) 

J2(B\b{a*(l-^)+(3*(l +1 5 )}u\X) 
x 

(X\u^{L u (l - 7 5 ) + R u (l + 7 5 )} M |0)(2tt) 4 5\p B - p K - p x ) 
In the parton model approximation this is written as [using Eq. (39)] 
MM = -if K (K~\sY(l - l 5 )b\B) 

°(B\b{a*(l - 7 5 ) + + 7 5 )}(^ u + m u ) 



(43) 



7 M {L u (l-7 5 ) + J R u (l+7 5 )}«|0) 



d 3 p u 



(27r) 5 (p b -p K - Pu) 



(44) 
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Kinematically this term looks like the two body decay b — > Ku. 
Using the definition 

(B-|67"(l-7 5 M0> = zy B p£ ( 45 ) 
and the quark equation of motion, we have 

2 

(B~|5(l + 7 5 H0> = -»/b _* (46) 

m b + m„ 

and finally we can write (dropping the u quark phase space factor and the delta function) 

M X M\ + M\M 2 = 2f B f K [ReC{g +PB ■ p u + g^p K ■ p u } + ReD{g + m 2 B + g_p B ■ p K }} (47) 

9+ = U + /- (48) 

9- = U- f- (49) 
2m 2 

C = J^( a *L u -(3*R u ) (50) 

m b + m u 

D = [P*L u -a*R u )]2m u (51) 

For decay to K*~ one can write a similar expression 

M X M\ + M\M 2 = ARe{FL* u L u )m K *g K *f B + 

i 

4Re(FL* u R u )m K *g K *f B (52) 

where the expressions for Xi, yi are given in the Appendix . 

There will be higher order perturbative corrections, such as additional gluons in the final states, 
to the processes described above. These effects are expected to be small. We will, however, include 
the bound state effects of the b quark inside the B meson on the decay distributions . 

Inside the B meson, the b quark is not on-shell. This will cause the energy to have a distribution 
even in the case of "two body" decay. To obtain the decay distribution we consider a model for 
the i?-meson wavefunction pH| which has been used for the calculation of the photon spectrum of 



inclusive b — > 57 decays in Ref. [27 



In the rest frame of the B-meson, the b-quark and the light antiquark q inside the B meson with 
energies E b and E q satisfy 



E b + E q = 09 s + m? q + ^p 2 + m 2 b (53) 
= m B (54) 



where p b = -p q , \p b \ =p = \p q \ 
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To satisfy Eq. (53) for all p the 6-quark mass is considered a function of p, 

m t = m B + m 5 — 2mBE q (55) 
The i?-meson wave function is taken as 

m = 7^ e_% (56) 

with the normalization 

roc 

p 2 (p(p)dp = 1. (57) 





For our numerical results we will use m u ^ = 150 MeV and the Fermi momentum pp = 0.3 GeV. 
The decay distribution is now obtained from 

,/J r p 2 dp<P( P )(^-) (58) 



cLEk J \ dEi< 1 

Ji \ xl / partonic 

where the partonic distribution (-jj^) partonic can be obtained by boosting the decay distribution in 
the rest frame of the b-quark to the rest frame of the B-meson. 

To complete the numerical calculations we have to fix the value of the gluon momentum q 2 in 
the G function of Eq. (4). For the "three body" decays governed by jM^ 2 , q 2 = p\ while for 
the "two body" decay governed by |M 2 | 2 and the interference terms one can use simple two body 
kinematics to obtain q 2 pa ml/2. In our calculation of M2 we compare results with q 2 = ml/ 3 
to those with q 2 = ml/2 in order to assess the dependence of the final result on this uncertainty. 

Having obtained the decay distributions we define the asymmetry for a B decay, 



^B -> K+^X) - £- k (B -> g-WjQ 



Integrated decay rates and integrated partial rate asymmetries can also be obtained in the usual 
manner. 

Following Ref [^] we can write the amplitudes, both Mi and M 2 as 

M t = V U A\ + V C A\ + V t A\ 

= V u Al + VAlt (60) 

where i — 1,2 and 

Kt = K-A\ 

Kt = K-a\ (61) 
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and the unitarity relation of the CKM, V u + V c + V t = has been used. Note that A l ut also contains 
the tree level amplitude. 

The decay distributions are such that the contribution from Mi and M 2 are separated in Ex 
in the approximation of neglecting the interference term. The decay distribution is a sum of two 
independent decay distributions governed by Mi and M 2 . From the structure of Mi and M 2 it can be 
seen that the Wilson coefficients that occur in Mi and M 2 contribute in pairs of the type Ci + c i+ i/N c 
and Ci + i + Ci/N c , respectively. The values of the Wilson's coefficients are such that generally the first 
combination is suppressed relative to the second and hence M 2 is enhanced relative to Mi. Thus 
the decay distribution associated with M 2 is larger than the decay distribution associated with M\. 
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FIG. 3. Predicted rate for B~ — > K~ X as a function of the kaon energy. Three curves in each figure 
are shown for N c = 2, 3, oo and provide an estimate of the theoretical uncertainty from the factorization 
hypothesis. In (a) and (b) different sets of form factors are then compared in order to determine the 
sensitivity of the predicted rate to the choice of a form factor model. The vertical scale in the plots is 
multiplied by (G F /V2) 2 x 1(T 6 . 



Using the form of Mi and M 2 given above one can write the partial rate asymmetries as 



2Im(V u *V c )Im(A i ut Ai t 



Eij (|K| 2 1 Ai t Ai t * | + |K| 2 |A^A J ct *| +2Re(V u *V c )Re(Ai t Ai t ' 



where i,j = 1, 2. The net asymmetry, a, is given by the sum 



a n + a 22 



.12 



,21 



(62) 



(63) 
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From the values of the Wilson coefficients c\ and c 2 it can be shown that the contribution to 
the asymmetry due to the interference of the tree and penguin amplitudes is suppressed in a 11 . 
This coupled with the fact that the gluon momentum q 2 is varying for M\ while it is more or less 
fixed for M 2 can lead to a larger value for a 22 compared to a 11 and a 12 . It should be pointed out 
that the interference term between M\ and M 2 can be important when calculating the partial rate 
asymmetries. 



VII. RESULTS AND DISCUSSION OF THEORETICAL UNCERTAINTIES 



In this section we discuss the results of our calculations which are shown graphically in Figs. |3| 
i"5] . We find that there can be significant asymmetries in B — ► K(K*)X decays especially in the 
region Ek > 2 GeV which is also the region where an experimental signal for such decays can be 
isolated. The branching ratios are of order O(10~ 4 ) which are within reach for future B factories. 



p=-0.15 

T)=0.23 

H=0.33 

ti=0.43 




p=0.15 



T|=0.23 
T]=0.33 
TH>.43 




2.1 2.2 2.4 2.5 2.7 2.9 3.0 2.1 2.2 2.4 2.5 2.7 2.9 3.0 



E K (GeV) E K (GeV) 

FIG. 4. Sensitivity of the predicted rate for B~ — > K~ X to p, rj. The three curves indicate the 
sensitivity of the rate to the value of Wolfenstein parameter r\ for fixed values of p 

The contribution of the amplitude with the top quark in the loop accounts for 60-75% of the 
inclusive branching fraction. However, since the top quark amplitude is large and has no absorptive 
part in contrast to the c quark amplitude, the top quark contribution reduces the net CP asymmetry 
from 30-50% to about 10%. This calculation includes the contribution from electroweak penguins. 
We find that the electroweak penguin contributions increase the decay rates by 10-20% but reduce 
the overall asymmetry by 20-30%. 

Let us now identify the main sources of uncertainties in our calculation. These are the use of 
the factorization approximation, the choice of a form factor model, the choice of q 2 for the gluon 
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FIG. 5. Sensitivity of the asymmetry for B — > K X to the Wolfenstein parameter rj for two fixed 
values of p. The three curves indicate the sensitivity of the asymmetry as a function of kaon energy. 



momentum in the G function in Eq. (4) for the "internal" two body diagrams, and the choice 
of a model for the 5-meson. We now discuss the sensitivity of the results for decay rates and 
asymmetries to these theoretical uncertainties. 

We have used the factorization approximation. The factorization approximation is expected to 
be valid in our calculations as we have argued at the beginning of Section 2. To take into account 
corrections to this approximation we allowed the number of colors to be a free parameter. In our 
calculation we consider the cases N c = 2, 3, oo although the analysis of exclusive two body B 
decays suggests that iV c ~ 2. In Fig. [3] we show the decay distribution for B~ — > K~X where X 
does not contain any strange particles. In the region of interest to experiment (i.e. Ex > 2.0 GeV) 
the decay distribution has only a modest dependence on N c . 

The second source of uncertainty is the choice of form factors used to describe the B — > K(K*) 
transitions. As mentioned earlier we use the form factors given in Ref . |23| . In Ref [f23] a monopole 
k 2 dependence for the form factors is chosen and the form factor at k 2 = for the D decays is 
fixed from semileptonic D — > K decays. The form factors are then scaled to the case of B — > K 
decay using heavy quark effective theory (HQET). The primary effect of choosing a different set of 
form factors in our calculations is to change the decay distribution of the "three body distribution" 
which is governed by the matrix element M\. Figures. 3(a), 3(b) show the effect of choosing two 
sets of form factors from Ref [2<J and Ref |24| with a different k 2 behavior of the form factors. For 
the energy region Ek > 2 GeV the form factor effects are negligible for the decay distributions as 
well as the asymmetries. For the remainder of the discussion and in the rest of the figures we will 
only show the results with form factors from Deandrea et al. p3| . 
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FIG. 6. Predicted asymmetries for B~ — > K~X and B~ — » K*~X as a function of the kaon energy. 
The value of N c = 2 is fixed. The two sets of curves indicate the sensitivity of the asymmetry to the values 
of q 2 for the gluon in the internal two body diagram. 



The asymmetries are sensitive to the values of the Wolfenstein parameters p and 77. The existing 
constraints on the values of p and 77 come from measurements of | V^t, | / 1 V^^ | , ex in the K system 
and AMs d . (See Ref. [29] for a recent review). In our calculation we will use fs = 170 MeV 
and choose (p = —0.15,77 = 0.33). To determine the dependence of our results on 77 we will also 
consider three sets of representative values, (p = —0.15,7] = 0.23), (p = —0.15,7/ = 0.33) and 
(p = —0.15,7/ = 0.43). We will also consider the set of 77 values with p = 0.15. The dependence 
of the rates on p, 77 is shown in Fig. ^, while the sensitivity of the asymmetry to these parameters 
is shown in Fig. |5|. For fixed p, the asymmetry increases monotonically as 77 increases. The results 
suggest that measurement of asymmetries in inclusive decays will give useful information on 77 once 
the size of the theoretical uncertainties is reduced. 

There are several other sources of uncertainty. These are: (1) the choice of q 2 for the gluon 
momentum in the G function in Eq. (4) for the "internal" two body diagrams and (2) the choice of 
a model for the wavefunction of the 5-meson. The q 2 variation causes a small change in the decay 
distribution but a fairly significant change in the asymmetries. This uncertainty is illustrated by 
comparison of the two curves in Figs. |], [15|. 

The choice of the value of the charm quark mass m c , is also a source of uncertainty for the 
asymmetry. We have taken m c = 1.35 GeV for this calculation. Increasing the charm quark mass 
to m c = 1.6 GeV does not significantly modify the decay rates but reduces the asymmetry by 
about 30%. Since the Wilson coefficients are calculated to next to leading order, they have little 
sensitivity to the renormalization scale. However, the G function which enters into the calculation 
has a stronger dependence on renormalization scale. Varying p? from ml/2 to 2m\ changes both 
the asymmetry and the decay rates by ±10%. A lower renormalization scale corresponds to a larger 
asymmetry. The model of the B meson wavefunction, used to take into account the Fermi motion 
of the quarks, is yet another source of uncertainty. We have used the model of Ali et al. which has 
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FIG. 7. Predicted asymmetries for B~ — » K~X and B~ — > K*~X as a function of the kaon energy. The 
three sets of curves indicate the sensitivity of the asymmetry to the value of N c . The values N c = 2, 3, oo 
are considered. 



been previously used for the B — > X s 7 case p7| . Different models give somewhat different results 
for the decay distributions while the asymmetries are insensitive to the choice of model for the B 
meson wavefunction. 
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In Fig. § we show the decay distribution for K* in the final state while in Figs 
the variation of the decay distribution and asymmetry with different sets of p, rj for 7V C 
Fig. 11 we show the asymmetry for several values of N c using the form factors from Ref p3j and 



13| we show 
In 



2. 



ml/2.0. A similar variation of the asymmetry with q 2 as calculated for the K in the final state 
is also observed in this case as shown in Fig. 10. 

Turning to B° decays, only the Mi part of the matrix element contributes. In Fig. |TT] we show 
the decay distribution for various N c values using the form factors from Ref and q 2 = m^/2.0. 
In Fig. [12| we show the decay distributions of B° — > K~ for representative values of [p,i]). In 
Figs. [13, [15| we show the asymmetries as we vary (p, if) and q 2 . The variation of the asymmetries 
with N c is negligible in this case. 

In Table. | we give the branching fractions and the integrated asymmetries for the inclusive 
decays for different N, q 2 = ml/2, fs = 170 MeV, p = — 0.15,7/ = 0.33 and the form factors from 
Ref. |E3|. For the charged B decays we also show the decay rates and asymmetries for Ek > 2 (2.1) 
GeV as that is the region of the signal. 

The above figures show that there can be significant asymmetries in B — >• K^X decays, espe- 
cially in the region Ek > 2 GeV which is the region of experimental sensitivity for such decays. As 
already mentioned, our calculation is not free of theoretical uncertainties. Two strong assumptions 
used in our calculation are the use of quark level strong phases for the FSI phases at the hadronic 
level and the choice of the value of the gluon momentum q 2 in the two body decays. Other uncer- 
tainties from the use of different heavy to light form factors, the use of factorization, the model of 
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FIG. 8. Predicted rate for B~ — > K*~X as a function of the kaon energy. The three curves indicate 
the sensitivity of the rate to the value of N c . The values N c = 2, 3, oo are considered. 



the B meson wavefunction, the value of the charm quark mass and the choice of the renormalization 
scale \x have smaller effects on the asymmetries. 

The use of quark level strong phases at the hadronic level neglects the possibility of soft FSI. 
This neglect may be a better approximation in our inclusive case as opposed to exclusive modes. 



VIII. COMPARISON WITH PREVIOUS CALCULATIONS 



Inclusive direct CP violating asymmetries were calculated by Gerard and Hou . They obtained 
asymmetries much smaller than the results given here. The differences between these calculations 
and our results can be understood as follows. In the Gerard-Hou calculation only the three body 
quark level process corresponding to b — > sqq (i.e. M\ in equation (1)) was considered and the limit 
N c — ■> oo was taken. With these conditions we obtain an asymmetry of —1.7% which then agrees 
with their result. We obtain large asymmetries only in the kinematic regime dominated by the two- 
body process b — > K u. These asymmetries agree qualitatively with the asymmetry for exclusive 
modes such as B~~ — ► K~n° found in several recent calculations in both sign and magnitude. In the 
kinematic region of B~ — > K~tt° we find an asymmetry of —(8 — 14)% to be compared with and 
— (2 — 8)% found by Kramer, Palmer, and Simma |llj and —(3 — 9)% found by Kamal and Luo [11 
for the exclusive mode B~ — > K~tt . 

We have also verified that the asymmetry in the process b — ► sec has the opposite sign to the 
asymmetry in the sum of b — > suu, b — > sdd, b — > sss (as does the asymmetry in B — > DD) as 
required by CPT to ensure the cancellation between all modes. 
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FIG. 9. Predicted rate for B~ — ► as a function of the kaon energy. The three curves indicate 

the sensitivity of the rate to the values of the Wolfenstein parameters p, rj. Only the signal region is shown 
in the figures. 

IX. CONCLUSION 

We find significant direct CP violation in the inclusive decay B — > K~X and B — > K*~X for 
2.7 > E K (*) > 2.0 GeV. The branching fractions are in the 10~ 4 range and the CP asymmetries 
may be sizeable. These asymmetries should be observable at future B factories and could be used 
to rule out the superweak class of models. 
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FIG. 10. Predicted asymmetry for B~ — > K*~X as a function of the kaon energy. The three curves 
indicate the sensitivity of the asymmetry to the values of the Wolfenstein parameters p, rf. 




FIG. 11. Predicted rate for B° — > K X and B° — > X as a function of the kaon energy. The three 
sets of curves indicate the sensitivity of the rate to the value of N c . The values N c = 2, 3, oo are considered. 
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FIG. 12. Predicted rate for B° K X as a function of the kaon energy. The three curves indicate 
the sensitivity of the rate to the values of the Wolfenstein parameters p, rj. 
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FIG. 13. Predicted rate for B° -»• K*~ X function of the kaon energy. The three curves indicate 
the sensitivity of the rate to the values of the Wolfenstein parameters p, rj. 
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FIG. 14. Predicted asymmetries for B° -► K~ X and B° -► K*~ X function of the kaon energy. The 
three curves in each figure indicate the sensitivity of the rate to the values of the Wolfenstein parameters 
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FIG. 15. Predicted asymmetries for B° -► K~X and B° K*~X as a function of the kaon energy. 
The two sets of curves indicate the sensitivity of the asymmetry to the values of q 2 for the gluon in the 
internal two body diagram. 
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TABLE I. Integrated decay rates and asymmetries for B — ► K^X Decay 



Process 



Branching ratio (1.65 x 1CT 



Integrated asymmetry 



B~ 



B~ -► K~X 
♦ K-X(E K > 2.1GeV) 
K-X (E x > 2.1GeV) 



1.02,0.79,1.20 
0.81,0.74,0.77 
0.6,0.7,0.8 



-0.10,-0.11,-0.050 
-0.12,-0.12,-0.07 
-0.12,-0.12,-0.13 



B~ 
B° 



> K*~X{E K , > 2.1GeV) 
> 2.1GeV) 



1.37,1.24,2.30 
1.05,1.16,1.67 
1.05,1.16,1.39 



-0.11,-0.14,-0.11 
-0.14,-0.15,-0.14 
-0.15,-0.15,-0.16 
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APPENDIX A: 



For B~ -> K~X decay 



A\ V H^ = \g + \ 2 D(p B ,p B ) + \g^\ 2 D(p K ,p K ) + 2Re(g+g*_)D(p B ,p K ) 
A\ V H^ = 



where 



g + = /+ + /- 
g- = f+- f- 



and we have defined 



D(A, B) = A- Pl B -p 2 + B ■p 1 A-p 2 - Pl -p 2 A- B 

g^H^ = \g+\ 2 m 2 B + \g-\ 2 m 2 K + 2Re(g + g*_)p B ■ p K 

For B~ — > decays (with p = p B + Pk*, <1 = Pb — Pk*), 

H^A\ u = J2T t 



while 



T i = b i{- 2 Pi ■ P2{rn 2 B m 2 K * - (p B ■ p K *f) 

+2p B ■ Pk*D(pb,Pk*) - m K * 2 D(p B ,p B ) - m B D(p K *,p K *)} 



T 2 = b\ 



2pi -p 2 H h—D(p K *,p K *) 



m 



K* 



T 3 = b 2 3 FD(p,q) 
T 4 = b 2 FD(p,q) 

T 5 = 2b 2 b 3 [D(q,p) - xD(p K *,p)} 
T 6 = 2b 2 b 4 [-D(q,q) +xD(q,p K *)\ 
T 7 = -2Fb 3 b 4 D(p,q) 



H^A 2 ^ = Ac x c 2 {pi ■ p B p 2 ■ p K * - Pi • PbPi • Pk*} 



and 

g, u H^ = J2S l 
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51 = -2b x [m K ,m B - (p B ■ p K *) } 

5 2 = -3b 2 2 

53 = Fb\p 2 
S, = big 2 

5 5 = 2b 2 b 3 [q ■ p - xp K * -p] 

5 6 = 2b 2 b i [-q 2 + xp ■ q] 

5 7 = -2Fb 3 b 4 p ■ q 



(All) 



where 



_ q 2 + (Pk* ■ qf 



vtik* 



x 



Pk* ■ q 



(A12) 



Finally for K* decay the interference term is, with pi = Pb, P2 = Pk* and q = Pi — P2, 

xi = 2b 1 [m 2 B p u -p 2 -pi- PuP2 ■ Pi] 



x 2 = -b 2 
x 3 = -63 

X4 = —64 
yi = -3b 2 



2pi ■ p u + 



D(P2,P2) 



M 2 



D(pi+p 2 ,q) - j^-D(pi+p 2 ,p 2 ) 



-D(q, q ) + ^D( q ,p 2 ) 



V2 = h 



(P1+P2) ■ q - (pi +P2) -P2 



P2 ■ q 



M 2 



2 P2 ■ q 

q - q - p2 ^r 



V3 = h 

D(A, B) — A ■ pi B ■ p u + A ■ p u B ■ p B — p u ■ pi A ■ B . 



(A13) 
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